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ABSTRACT 



Aims. We investigate the signature of the photospheric albedo contribution in solar flare hard X-ray spectra, the effect of low energy 
cutoffs in electron spectra, and the directivity of hard X-ray emission. 

Methods. Using Ramaty High Energy Solar Spectroscopic Imager (RHESSI) flare data we perform a statistical analysis of spatially 
integrated spectra and positions of solar flares. 

Results. We demonstrate clear centre-to-limb variation of photon spectral indices in the 15-20 keV energy range and a weaker depen- 
dency in the 20 - 50 keV range which is consistent with photospheric albedo as the cause. The results also suggest that low-energy cut- 
offs sometimes inferred in mean electron spectra are an artefact of albedo. We also derive the anisotropy (ratio of downward/observer 
directed photons) of hard X-ray emission in the 15-20 keV range for various heliocentric angles. 

Key words. Sun: flares, Sun: X-rays, Sun: particle emission, Scattering, Techniques: spectroscopic 



1. Introduction 

The successful operation of t he Ramaty Hig h Energy Solar 
Spectroscopic Imager RHESSI dLin et al.ll2002l) since February 
2002 has provided us with a substantial database of solar hard 
X-ray flares that can be used for detailed spectroscopic analysis. 
The high accuracy of the spectroscopic data, in combination with 
the positions of solar flares, allows us to determine the centre-to- 
limb spectral variations in various energy ranges. 

Variation of hard X-ray flare occurrence and spectral in- 
dex with solar position of the sources have be en of great in- 
terest fo r severa l decad e s. As suggested e. g. by lElwert & Haug 
197 lb; iBrownl d 19721) : IPetrosianl (119731) : iLeach & Petrosianl 
1983), centre-to-limb variations could provide essential infor- 
mation on directivity of X-ray emission and hence acceler- 
ated electrons. Most earlier analyses concentrated on centre- 
to-lim b variation of flare occurrence. lOhkil d 19691) and iPinterl 
(1969) reported significant centre-to-limb variation at energies 
above 10 keV although they did not correct the hard X-ray 
flare o ccurrence for Hg flare longitude distribution. Subsequent 
works dDrakdll97UlCatalano & van Allenlll973l: lPhilUpslll973t 
Datlo we et al.ll974tlKanell974t) covering various energy ranges 
from 1-10 ke V claimed no centre-to -limb variation of hard X- 
ray occurrence. iDatlowe et al.l (1 19771) concluded that anisotropy 
in the 10 - 100 keV range was so low that it ruled out the 
downward-beamed thick targ et model du e to lack of limb 
brigh tening predicted e.g. by (Brown 1972). Th eoretical mod- 
els (IPetrosianl 119731: lLanger & Petrosianl [19771: iBai & Ramatvl 
1 1978b predicted that at hard X-ray energies, above 15 keV, 
centre-to-limb variations of spectral properties should be sen- 
sitive to elect ron directivity. Altho ugh such studies were not 
very common, IDatlowe et al.1 dl974 and lRov & Datlowd d 1975b 



Send offprint requests to: Jana Kasparova 



found, contrary to above mentioned works based on hard X- 
ray occurrence, significant centre-to-limb spectral variation in 
the 10 - 100 ke V range and attributed it to anisotropy of X-ray 
emission. Then, IVestrand et alJ (1987) demonstrated that Solar 
Maximum Mission (SMM) GRS data in the 25 - 200 keV range 
and above 300 keV, show significant spectral ha rdening of events 
near th e limb. Similar results were reported by lBogovalov et all 
d!985l) from Venera 13 data in the 50 - 100 keV range and by 
McT iernan & Petrosianl dl99ll) in the 0.3 - 1 MeV range. 

It has been known for a while that X-ray photons can be 
effectively backscattered by photo sphere atoms and electrons 
(TombHrJl97l:Bai & Ramaty |j978j). Thus, at energies not dom- 
inated by absorption the backscattered albedo flux must be seen 
virtually in every solar flare spectrum, the degree of the albedo 
contr ibution depending on the directivity of the primary X-ray 
flux dKontar et al.l l2006). The solar flare photons backscattered 
by the solar photosphere can contribute significantly (the re- 
flected flux is 50-90 % of the primary in the 30 - 50 ke V range 
for isotropic sources) to the total observed photon spectrum. For 
the simple case of a power-law-like primary solar flare spectrum 
(without albedo), the photons reflected by the photosphere pro- 
duce a broad 'hump' component. Photospheric albedo makes the 
observed spectrum flatter below ~ 35 keV and slightly steeper 
above, in comparison with the primary spectrum. The amount of 
backscattered photons is given by the scattering cross-section 
and is roughly dependent on the projected area of the albedo 
patch, which makes the photospheric albedo stronger for so- 
lar centre events. Photospheric albedo has been treated using 
Monte Carlo simulations to calculate the o bserved spectrum fo r 
an assumed primary power-law spectrum (Bai & Ramaty 19781) . 
iKontar et al.l ((2006) developed a new a pproach to photospheric 
albedo based on the Green's functions o fMagd ziarz & Zdziarski 
d 19951) that produces the correction for any solar flare X-ray 
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Fig. 1. Photon spectrum 1(e) (upper panel) and spectrum in- 
dex y(e) (lower panel) for a modelled primary photon spectrum 
I p (e) ~ e~ 3 for an isotropic source located at n — 0.8. Primary 
l p (e) and y p (e) are plotted with a dashed line and the observed 
I (e) and y a (e) with a solid line. Horizontal lines indicate spec- 
tral indices yo, y\, and 72 obtained by single power-law fits in 
the 15 - 20, 20 - 35, and 35 - 50 keV energy ranges, respectively. 



spectrum and can be used both for forward and inverse spectrum 
calculations. The influence of albedo was not taken into account 
for spectral centre-t o-limb variation in ear lier analyses and it was 
wrongly dismissed dDatlowe et al. 1974) a s being unimpo rtant 
in the 20 - 50 keV r ange. However. iBai & Ramatvl d 19781) and 
IVestrand et al.l dl987l) pointed out that the Datlowe's results are 
consistent with an al bedo contribution to t he observed photon 
spectrum. Using the IBai & Ramatvl dl978) results, the impor- 
tan ce of albedo i n electr on spectrum infere nce was first included 
Johns & Linl d 19921) and analysed by lAlexander & Brownl 



t hat the 
dBrownl 

197lt Brown et IDE003) can exhibit low-energy cutoff - e.g. 



On the other hand, it has been long believed 
mean electron flux distribution F(E) defined by 



Nitta et all d 19901) showed that photon spectra of the impulsive 



component flatten toward low energies and suggested that a low- 
energy cuto ff of the electron s pectrum around 50 keV could be 
responsible. iFarnik et alj d 19971) . using Yohkoh data, observed a 
flare with a flat spectrum belo w 33 keV (photon spect ral index as 
low as 1.98). More recently, Kas parova et alj d2005l) reported a 
clear low-energy cutoff inferred from the observed photon spec- 
tra of the August 20, 2002 flare , neglecting albedo. Usi ng a 
regularised inversion technique, iKontar & Brownl d2006al) re- 
port three flares with a similar feature - an absence of electrons 
around 20 - 30 keV between possibly thermal and non-thermal 
component. 



Howev er, applying the mod el-independent albedo correction 
method of IKontar et all d2006l) to the August 20, 2002 flare, 
Kasparov a et al.l d2005l) have clearly demonstrated that such a 
high value of low-energy cutoff is removed by albedo correc- 
tion for an isotropic source. Equally, the albedo corrected pho- 
ton s pectra do not require a gap in the inverted ele ctron spec- 
trum (Kasparo va et al.l20 05; Kontar & Brown 2006a). Using the 



present survey, we can potentially find many flares with such low 
spectral indices that they might have clear low-energy cutoff ap- 
pearing as a gap between thermal and non-thermal distributions. 

The spectrum-independent photospheric albedo correction 
method is introduced in Section 2. Data selection and analysis 
are discussed in Section 3. In Section 4, our survey shows that 
the centre-to-limb variations of spectral indices are consistent 
with photospheric albedo. Moreover, we find limits on the X- 
ray spectral directivity for various heliospheric angles. Section 5 
confirms the prediction that the low-energy cutoff in F(E) is re- 
moved when the spectra are albedo corrected for all events with 
low spectral index found in the survey. Lastly, Section 6 reviews 
the results obtained and their consequences for the angular dis- 
tribution of X-ray emitting electrons. 

2. Photospheric albedo and spectral index 

The albedo contribution produces substantial change to the ob- 
served spectrum in the range between 10 and 100 keV. It be- 
comes negligible for energies below 10 keV due to the photo- 
electric absorption by atoms while photons above 100 keV pen- 
etrate deep into high density layers of the solar atmosphere and 
do not escape back. Therefore, the reflectivity of the photosphere 
has a broad maximum in the 30 - 50 keV range. The total spec- 
trum in the observer direction [i(n = cos 9, where 9 is the helio- 
centric angle) can be expressed as 



I (e,fi) = I p (e,fx) + Ia(£,m) 
7 A (e,yu) = G(e, e' , fj)I p (e' , fx)a(ji) 



(1) 



where I p (e,fi) is the primary photon spectrum of the flare, 
7a(£, AO is the spectr um reflected into th e observer direction, G is 
the Green's matrix (IKontar et al. 2006) accounting for Compton 
backscattering and photoelectric absorption, and a(fj) is a pa- 
rameter of emission anisotropy, roughly the ratio of the flux to- 
ward the Sun I p (e,/j < 0) = a(p)I p (e, fi) and the flux toward the 
observer I p (e,fj). We simplify the problem here by consi dering 
only a(fj) independent of e (cf I Alexander & Brownll2006l) . 

The photon spectrum, lA(e,fi), backscattered toward the ob- 
server strongly depends on three parameters: the spectral index 
of the primary spectrum, the location of the flare on the solar 
disc, and the directivity of X-ray emission (anisotropy parame- 
ter a(p)). The contribution of the albedo to the observed spec- 
trum is greater for flatter primary spectra (low y) and for solar 
disc flares (large fi). Further, the more X-rays are beamed down- 
wards a(fj) > 1, the larger is the contribution of the albedo to the 
observed spectrum. 

The energy dependent spectral index 



7(e) = ~ 



e 61(e) _ dln(/(e)) 
1(e) de ~ dln(e) 



(2) 



for the primary spectrum y p (e) substantially differs from that 
(7o(e)) of the total I (e, fi) toward the observer (including albedo) 
- see Figure 1 for the case of y p = 3. In the range 10-25 keV 
the photon spectrum becomes flatter, y < y p , while the observed 
spectral index is higher at energies above 35 keV, y > y p . 
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Using a single power-law fit in various energy ranges 
(Figure [1), we can estimate the local spectral index assuming 
it is constant in each energy range. (Note that generally for a 
photon spectrum 1(e ) ~ e~ c(e) , c(e) = 7(e) only for c = const 
dConwav et ai]|2003l) .) 

Hence, as a result of heliocentric angle dependent albedo, 
the spectral shape of photon spectra should vary as a function of 
their position on the solar disc, i.e. a centre-to-limb variation of 
spectral indices is expected and the degree of the spectral index 
variation with heliocentric angle could usefully provide us with 
the directivity of hard X-ray sources, i.e. the anisotropy coeffi- 
cient a(iS). 

3. Data selection and analysis 

In our analysis we used the whole RHESSI flare list (from 
February 12, 2002 until April 18, 2006) to search for any evi- 
dence of centre-to-limb variation in photon spectra. 

RHESSI has 9 detectors on board with FWHM resolution 
about ~ 1 keV. For X-ray spectroscopy analysis we use the 
six front segments of the detectors ignoring detector 2, 5 and 
7 due to either l ow resolution or the segmentation problem 
(ISmith et al]l2002l 

To handle such a large set of data, we set up an automatic 
procedure which selects suitable events and retrieves the flare 
photon spectra. To be chosen for further analysis the events had 
to meet the following major criteria: 

- The spectra of the events should be observed above 25 keV. 
Thus, we have selected flares flagged as observed above that 
energy from the RHESSI flare list and have determined the 
time of peak flux in the 25 - 50 keV range. The events flagged 
as particle events were ignored. This led us to ~ 1500 events. 

- Count rates exceeding about2xl0 3 counts per second per de- 
tector cause pulse pi le-up in the detect ors that requires trou- 
blesome correction (Smith et al. 2002). To avoid the pile-up 
issues, we have thrown away events with corrected livetime 
counter < 90%. We disregarded also weak events (back- 
ground subtracted count rate lower than 5 counts per second 
in the 25 - 50 keV range) and events possibly affected by X- 
ray absorption in Earth's atmosphere (peak time closer than 
60 s to RHESSI eclipse time). This reduces our survey to 
~ 800 events. 

- We have used CLEANed images in the range 10 - 20 and 
20 - 50 keV to check for true solar events and to find their po- 
sitions. Finally, the total number of the selected flare events 
drops to 703. 

RHESSI is an instrument with high background that makes 
background subtrac tion an important step in the data analysis 
dSchwartz et alj[2002b . The code searches for two time intervals, 
one before and one after the flare peak, where count rates are 
less than 0.01 of the peak count rate in 15 - 50 keV range and 
0. 1 above 50 ke V. The background time intervals are determined 
for several energy bands separate ly and the backgro und at the 
peak is estimated by a linear fit dSchwartz et al.ll2002l) . See the 
flare example in Figure[2] 

The count spectra [counts sec -1 keV 1 cirT 2 ] were accumu- 
lated in 12 s time intervals with a quasi-logarithmic energy bin- 
ning from 10 to 100 keV to improve the signal-to-noise ratio 
(SNR) at higher energies (Figure [2j. After background sub- 
traction, the resulting peak count flux (one spectrum per flare) 
and the detector response matrix (DRM) for the instrument at 
this time inte rval were used as an input for t he regularised inver - 
sion routine dScullion et alJIin preparation! iKontar et"aT1 12004). 
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Fig. 2. Example of the flare spectrum for July 30, 2002 flare, 
17:37:36-17:37:48 UT. The upper panel corresponds to the reg- 
ularised photon spectrum, the lower panel shows the count spec- 
trum. Observed spectrum is plotted with thin crosses, back- 
ground subtracted flare spectrum with thick crosses. The back- 
ground is displayed as a histogram. Horizontal sizes of crosses 
and width of the histogram bins correspond to widths of energy 
bands, while vertical sizes represent lcr statistical uncertainties 
of the flux. 



The regularised inversion of a count flux spectrum gave us a non- 
parametric photon spectrum and its uncertainties. 

3. 1 . Photon spectra of selected events 

As discussed in Section [2] addition of the albedo contribution 
changes the shape of the primary spectra. To describe its ef- 
fects, the regularised photon spectra have been fitted with a sin- 
gle power-law in three separate energy ranges 15 - 20, 20 - 35, 
and 35 - 50 keV thus obtaining corresponding spectral indices 
7o,7i, and 72, respectively. 

Note that spectra in the 15-20 keV energy range may be 
affected by thermal component(s) which could be compara- 
ble to or dominate the non-thermal emission and masking 
the albedo contribution. Considering a widely used isother- 
mal approximation, value of 70 in such flares would be higher 
in comparison with events without strong thermal compo- 
nent. 

Since the count flux decreases with energy while the back- 
ground level does not (Figure |2), naturally not all 703 events 
have significant flux above background at all three energy 
ranges. For 70 an d J\ analysis we accepted only those events 
with count flux signal-to-noise ratio (SNR) > 3 at all three en- 
ergy ranges. Furthermore, the Green's matrices and therefore the 
corrected DRM are non-diagonal at 72 energies, the photon flux 
at this energy range depending on the signal at higher energies 
and making 72 unreliable for flares with low SNR above 50 keV. 
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Thus, for j2 analysis additional condition SNR > 3 in the 50 - 
66 keV energy range was required. 

Hence, the above conditions select only events with signif- 
icant flux above the roughly estimated background in energies 
higher than the energy range where the spectral indices are de- 
termined. This reduced our set to 398 events for yo and j\ analy- 
sis and to 123 events for j2 analysis - see Figure[3] Note that the 
events are not uniformly distributed over fi. The main reason is 
that active regions and flares are located in a rather small range 
of latitudes. Thus fewer events per unit fi are observed close to 
the limb than at the disc centre. This also puts a constraint on 
the width of bins for anisotropy analysis - see Section l4~T1 and 
Tabled 

We also apply an isotropic (a(fi) = 1) albe do correction 
using Equation ([T]i and the approach discussed in iKontar et al.l 
(2006) to determine spectral indices of albedo corrected spectra. 
The relative occurrence of spectral indices of observed spectra 
and of those corrected for isotropic albedo is shown in Figure [3] 
Note the excess of events with yo < 2 which is removed after the 
isotropic albedo correction. The excess of yo > 7 corresponds to 
events with very steep spectra. 
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4. Centre-to-limb variation 

The resulting spectral indices yo, y\, and 72 of the flares are 
shown Figure [4] There is no clear upper bound on y val- 
ues (primary and observed spectral index can be arbitrar- 
ily large) and the albedo contribution to the observed spec- 
trum becomes negligible for steep primary spectrum - see 
Section |2j On the contrary, the primary spectral index for 
bremsstrahlung emission cannot be lower than 1. Therefore, 
for display purposes only values up to 3 (3.5 for 72) are shown; 
events with larger spectral indices fill the whole range of fi. 

The spectral index yo shows significant centre-to-limb varia- 
tion (see Figure 0). Lower values of yo tend to be located closer 
to disc centre, forming an edge in the yoip) distribution. There 
are no flares with yo < 2 located at yU < 0.5 (9 > 60°), all such 
being at > 0.5. We do not see similar pattern at higher energies 
(spectral indices y\ and 72) and the Figure [4] shows no correla- 
tion with the position of flares in these energy ranges. 

As one can see from our albedo model with a + a(e) 
(Figure Q}, the strongest centre-to-limb variation due to photo- 
spheric albedo is expected at lower energies, i. e. for yo, and al- 
most no dependency is expected on flare position at higher ener- 
gies. The corresponding theoretical model for single power-law 
primary spectra is shown by various lines in Figure|4] 

Figure [4] also shows that anisotropy plays a significant role 
in centre-to-limb spectral variation only at yo energies. Although 
we have no information about the primary spectra and distribu- 
tion of primary spectral indices of analysed events, it is clear 
from the figure that the albedo model predictions, assuming sin- 
gle power-law primary spectra, are consistent with the general 
variations of the observed spectral indices with /i. In order to put 
further constraint on the model, particularly at energies above 
20 keV, spectral indices need to be determined with uncertain- 
ties smaller than the predicted centre-to-limb variation - see 
Figure |4] 

Following prev i ous results ( Ves trand et al.l 119871; 
iDatlowe et all 11974 1 19771) . and to assess the change of 
spectral indices with jj, we divide the sample into two subsets; 
jj < 0.5, and /j. > 0.5 and for each interval we calculate the mean 
values of yo, yi, and y2 of observed and isotropic albedo cor- 
rected phot on spectra. We the n apply the Kolmogorov-Smirnov 
(K-S) test (iPress et all 1 19921) to calculate the probability that 
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Fig. 3. Top: Positions of selected flares. Crosses indicate 398 
flares used for the centre-to-limb analysis of yo and y\, small 
circles correspond to 123 flares for the y2 analysis (see text for 
details). Large circles indicate regions of /i = 0.2, 0.4, 0.6, and 
0.8. Middle and bottom: Relative occurrence of yo (solid), y\ 
(dashed), and y2 (dotted) from observed spectra and those cor- 
rected for isotropic albedo, respectively. 
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the distributions of spectral indices at fi < 0.5 and fi > 0.5 
are drawn from the same parent distribution. The mean values 
and the probabilities (see Table Q]) show the trend expected for 
isotropic albedo. The mean yo value of yo increases after the 
albedo correction is applied as one would expect. The mean 
spectral indices y~i, y~2 with and without albedo correction do 
not differ significantly - i.e. are within the uncertainties of the 
mean values. The K-S probabilities show a similar trend, the 
distributions with the albedo component removed show larger 
probabilities of being drawn from the same dataset. However, 
the probability values for observed indices are too large to reject 
the null hypothesis, partly because the albedo contribution is 
smeared over too broad a range of fi, Afi = 0.5. 

Better results can be achieved by taking two more extreme 
distinct sets of flares; one close to the limb with fi < 0.1 (6 > 
84°), where the albedo is negligible, and disc centre events with 
fi > 0.9 (9 < 25°) where the albedo should be strongest. The 
value of K-S probability for the yo distributions drops to 0.9%, 
allowing us to reject the null hypothesis at 0.05 and 0.01 sig- 
nificant levels and conclude that the yo distributions are 'sig- 
nificantly different'. 

After applying the isotropic albedo correction, the yo dis- 
tributions can no longer be considered as being different - see 
Table Q] and the yo distributions in Figure [5] However, the same 
conclusion cannot be reached for yi and j2 - see Table Q] It 
should be pointed out that the number of events for j2 analysis 
for those fi bins is rather low - see Table Q]- so a larger set of 
events will be needed to reach any conclusion. 

4.1. Directivity of X-ray emission 

The intensity of photons backscattered from the photosphere is 
determined by the downward directed primary photon spectrum. 
The larger the downward directivity of the source, the stronger 
is the albedo contribution to the observed spectrum (see the pre- 
dictions for a = 1,2,4 and the model of a single power-law 
primary spectrum in Figure U)- We describe the downward ver- 
sus toward observer flux ratio by the anisotropy parameter a(fi), 
see EquationQ] (which we take to be independent of energy). 

In order to assess the range of primary X-ray directivity con- 
sistent with the data, we have assumed limb spectra, I {e,fi < 
0.1), to be true primary spectra toward the observer, I p (e,fi) = 
I (e,fi < 0.1), and added to these the albedo component for var- 
ious a(fi). Using Equation (Q]l, the modelled spectra Im in the 
observer direction fi take the form 



(3) 



7 M (e,/0 = [1 + a(n)G(£,e',fi)]I ( £ ',fi < 0.1) . 

This approach does not decrease the signal-to-noise ratio of 
photon flux as in the case of albedo subtraction which reduces 
the photon flux but keeps the corresponding uncertainty at the 
same level (uncertainty on the albedo contribution is set to zero.). 

The K-S test was then applied to distributions of spectral 
indices of the modelled Im and observed I a spectra in several 
ranges of fi. The width of the fi bins was chosen to contain ap- 
proximately the same number of events. Due to the non-uniform 
distribution of events with fi, see Figure[3] the width of the fi bins 
increases towards the limb (smaller fi). We regard the directivity 
a(ju) as consistent with the data if the K-S probability is larger 
than 0.05 or 0.01 significance level, i.e. when the null hypothesis 
that two distributions of observed and modelled spectral index 
for given a(fi) are drawn from the same distribution cannot be 
rejected. In the yo energy range we find the ranges of allowable 
directivity to be 0.5 < a(p) < 4 and 0.2 < a(ji) < 5, respectively 




0.5 



0.0 


0.2 


0.4 


0.6 


0.8 1 


3.0 






/j,=cos6 




I 1 

. + 


< 

O 


r IV i | 

O 

O 


i<7 i | i i 
<X> 

* oo 
o 


4- - 


2.5 


> - 


o 


O 





/>+ - 

- o+ *- 








o 


o 




2.0 
























1 .5 


1 . 


, , 1 


, , , 1 , 


, , 1 , , 





0.0 0.2 0.4 0.6 0.8 1.0 

/X = COS0 



2.5 - 




.0 - — 



1 .5 



0.0 0.2 



0.4 0.6 

^ = 0030 



0.8 



.0 



Fig. 4. Spectral indices yo, j\ , and j2 versus cosine fi of helio- 
centric angle 0. Vertical error bars indicate average uncertainties 
on the values as determined from single power-law fits. Lines 
show the predicted dependency for single power-law primary 
spectra with y p = 2.0,2.5 and for a(fi) = 1,2,4 (solid, dashed, 
and dotted lines respectively), see also Section l4~Tl High values 
of spectral indices are not shown to emphasize their lower 
bounds. Crosses represent flares for which F(E) was deter- 
mined from the total spectra. Flares with a dip in F(E) are 
denoted as stars (see SectionH). 
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4.2 ±0.1 (81) 


3.6 ±0.1 (29) 


K-S prob. 






0.4 


0.4 


0.8 



Table 1. Mean spectral indices yb, j\ and j2 for all selected events at selected ranges of p. K-S prob. corresponds to probability that 
a pair of distributions comes from the same parent distribution. The number of events in each p bin is indicated in brackets. 




Fig. 5. Distributions of yo for 'centre' (p > 0.9, dashed line) and 'limb'(/i < 0.1, solid line) flares. The left plot corresponds to 
observed yo, the right plot is for spectra corrected for isotropic albedo, i.e. with directivity a = 1. 



- see Table [2] and Figure [6] The results are also consistent the 
hypothesis that a(p) does not vary significantly with p for the 
0.05 and 0.01 significance levels. 

In the energy range of y\ no conclusions about anisotropy 
can be drawn because high values of the K-S probabilities for the 
distribution of y\ do not permit us to distinguish them, though, 
the results are not inconsistent with the albedo model which does 
not predict significant variation of j\ with p and a in this energy 
range - see Figure [4] 

Finally, the small number of suitable events at high energy 
prevented us from analysing j2- On the basis of the albedo 
model, see Figure[4] we do not expect to find a significant influ- 
ence of the anisotropy in the j2 energy range, but this prediction 
is not tested by the presented data. 
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5. Low-energy cutoff in F(E) demanded by the data 

The survey presented in the this paper is also suitable for analysis 
of flares where the hard X-ray spectrum may indicate suspicious 

features in F(E). 

Using a regularised inver sion technique dKontar et al.| [2004; 
Scullion et al. in preparation), we determined the mean electron 
flux distribution F(E) from count spectra of several events which 
had either flat photon spectra or a large change in local spec- 
tral index with energy. Such spectral behaviour is expected for 



Fig. 6. Directivity at different p determined from the distribu- 
tions of yo at the limb and a given range of p. Cross-hatched 
areas show the directivity for which the hypothesis that the ob- 
served distribution at a given range of p, and the distribution /m 
at the limb with added albedo contribution, are drawn from the 
same parent distribution cannot be rejected at the 0.05 signif- 
icance level; crossed areas correspond to the 0.01 significance 
level. See also Table|2] The dashed line shows the isotropic case, 
i.e. a(p) = 1. 
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K-S prob. for a(fi) 







N 


10 


5 


4 


2 


1 


0.5 


0.25 


0.2 


0.1 - 


0.4 


80 


3 x 10~ 3 


0.03 


0.07 


0.2 


0.07 


0.02 


8 x 10~ J 


6 x 10~ 3 


0.4- 


0.7 


88 


9 x 10~ 6 


4 x 10~ 4 


1 x 10- 3 


0.02 


0.4 


0.1 


0.03 


0.02 


0.7- 


0.9 


100 


7 x 10- 6 


2 x 10- 4 


4 x 10" 4 


0.02 


0.3 


0.04 


7 x 10- 3 


7 x 10~ 3 


0.9- 


1.0 


81 


2 x 10~ 3 


9 x 10~ 3 


0.04 


0.2 


0.5 


0.06 


0.01 


4 x 10~ 3 



Table 2. K-S prob. for several values of directivities at different ranges of /u (N is the number of events in that range). The values 
correspond to probability that the observed distribution of yo at a given range of fi and the modelled distribution Im of yo determined 
from the events at the limb (y < 0.1) with added albedo contribution for the given ju and directivity a are drawn from the same 
parent distribution. 



photon spectra either affected by the albedo - see Figure [TJ - 
or produced by F{E) with an absence of electrons in same en- 
ergy range. We found several new cases which exhibit a statis- 
tically significant dip in F(E) - e.g. Figure [7] (left). The posi- 
tions and spectral indices of events with such a dip are indicated 
in Figure |4] However, when the isotropic albedo correction was 
applied (i.e. a{p) = 1), all such dips in F(E) were removed - 
e.g. Figure [7] (right) . Our set of events includes the Sep 17, 2002 
(Kon taret all 120061) and Apr 25, 2002 flares dKontar & Brown! 
2006a]) for which a statistically significant dip in F(E) at ener- 
gies below 30 keV has been reported, neglecting albedo. 

The albedo contribution in events close to the limb should 
be negligible. Therefore, features in F(E) for events there could 
not be ascribed to the distortion of the photon spectra by albedo. 
F(E) was determined for a number of events close to the limb 
which showed flattening of the photon spectra. None of these 
events close to the solar limb (fi < 0.5) was found to have a 
statistically significant dip in F(E), but revealed only some flat- 
tening of F(E). All events with determined F(E) are indicated 
in Figure |4j 

6. Discussion and conclusions 

The high energy resolution of RHESSI down to ~ 1 keV and a 
large database have allowed us to study statistically 398 events 
not affected by pulse pile-up but having sufficient count rate to 
measure local spectral indices in three different energy bands. 
Combining these with flare position measurements we analysed 
centre-to-limb variation of hard X-ray flare spectra. The data 
clearly show a distinct variation of low energy spectral indices 
with heliocentric angle. This dependency is strong for spectral 
index in the range 15-20 keV but very weak for higher ener- 
gies. The mean values of low energy index yo are decreasing 
for larger y. (smaller heliocentric angles). We have compared 
centre-to-limb spectral variatio ns with the predictio ns of albedo 
induced spectral index change (Konta ret al.ll2Q06h for primary 
spectra assumed to be power-law and found the behaviour of 
mean spectral indexes yo, j\ , and yo (energy ranges 15-20 keV, 
20 - 35 keV, and 35 - 50 keV) to be consistent with the spectral 
change due to Compton backscattering albedo. For y\ and y2 at 
higher energies there is very little constraint information in the 
data. 

The intensity of backscattered photons is determined by in- 
tensity of downward beamed photons, so we can estimate the ra- 
tio of downward/observer directed photons. Although our data 
set contains large number of events of photon spectra with 
high energy resolution, it can be only used to test the hard X- 
ray anisotropy below 20 keV. Even for an anisotropic primary 
source, the albedo induced spectral change and centre-to-limb 
variation are significant only below 20 keV - see Figure QJ and [4] 



(which applies for a primary power-law photon spectrum). The 
inferred anisotropy ratio, 0.2 < a(p) < 5, (Figure|6]l in the 15 - 
20 keV range, and especially its large spread, allows no clear 
conclusion about the anisotropy. Moreover, at these energies 
anisotropy of electron beam emission can be masked by isotropic 
thermal emission, thus making the determination of anisotropy 
from spatially integrated spectra difficult. 

The value of anisotropy ratio obtained is consistent with 
the predictions of hard X-ray emission of downward colli- 
mated beam and cannot reject even qui te strongly beamed cases. 
For example, the detailed model of iLeach & Petrosianl d!983l 
Fig.4) predicts a 55 3 at 22 keV for disc centre events. On 
the other hand, the determined anisotropy a(jj) is also con- 
sistent with isotropic X-ray emission and lends support to 
quite isotropic electron d istribution when put together with the 
iKontar & Brown! d2006bi) analysis of detailed spectra of individ- 
ual events. The anisotropy parameter a(y) was also found not to 
need significant dependence on y. - see Figure [6]- indicating that 
the beam X-ray emission need not strongly vary in the upward 
hemisphere (different values of fi probe different viewing angles 
of upward emission). We stress that with a(fi) we measure the 
anisotropy of downward directed photons rather than upward di- 
rected photons at different angles as was done in most previous 
studies of solar flare directivities. 

The data analysed in the paper, as well as p r evious published 
results dNitta et al.l [l990l: iFarnfk et alj Il997t iKasparova et al.l 
2005; IKontar & Bro wn 2006a) describe events with observed 
photon spectra which are flat at low energies. Such behaviour 
suggests a low-energy cutoff in the mean electron spectrum as 
derived from such observed photon spectra. We show that all 
such flares were located close to the solar disc centre (jj > 0.5 
(9 < 60°)), where the albedo c omponent must be properly taken 
into account. [Nitta et ai1dl990l) reported that two flares observed 
with SMM on 15 October 1981 at 4:43 UT and on 29 March 
1980 at 9:18 UT had y ~ 2 and y = 1.9, respectivel y, both flares 
being located not far from disc centre (yu > 0.6). iFarnik et al.l 
d 19971) observed a few flares with very flat observed hard X-ray 
photon spectra from t he same active region near the disc centre 
(jj. = 1.0). We note fhat lZhang & Huangjd2003Ll2004 studied the 
combined effect of the low-energy cutoff and albedo contribution 
to explain very flat photon spectra observed by Yohkoh/HXT and 
concluded that both effect could cause flattening of photon spec- 
tra. However, to account for albedo they used integral reflectivity 
(angle-averaged) which is not applicable to centre-to-limb stud- 
ies. 

We have applied an isotropic albedo correction to the flares 
in our survey and have found this makes the alleged low-energy 
cutoff (or a 'dip') in electron spectra statistically insignificant. 
This survey strongly suggests the hypothesis that the low-energy 
cutoff inferred in the mean electron flux distribution is an artefact 
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19-JjI-2004 20:56:52.00 19-Jul-20C4 20:56:52.00 




10 100 10 100 

E [keV] E [kcV] 



Fig. 7. F(E) corresponding to the total (left) and primary (right) flare photon spectrum for an event close to the disc centre (19 Jul 
2004). Vertical error bars correspond to 3<x uncertainties, horizontal lines show the FWHM of energy resolution determined from 
regularised inversion. 



of the photospheric albe do. Hence, we ca n confidently confirm 
our earlier conclusion (iKontar et al.ll2006l) that low energy cut- 
off is a feature connected with albedo and not a physical prop- 
erty of the mean electron distribution. This result can substan- 
tially change the estimated flux and total number of electrons 
accelerated in a flare. Therefore, we strongly suggest that the 
albedo component must be considered in determination of elec- 
tron beam properties from spatially integrated solar hard X-ray 
spectra. 

Acknowledgements. We are grateful to H. S. Hudson and C. M. Johns-Krull 
for their valuable comments. JK acknowledges the financial support of a Royal 
Society Incoming Short Visit grant and a PPARC RA post to visit University 
of Glasgow in 2005 and 2006, and of grants 205/04/0358 and 205/06/P135 of 
the Grant Agency of the Czech Republic, and of research plan AV0Z 10030501 
of the Astronomical Institute AS CR. EPK and JCB acknowledge the financial 
support of a PPARC Advanced Fellowship and Rolling Grant respectively. This 
work was also supported by the Visiting Science Program at International Space 
Science Institute, Bern, Switzerland. 



References 

Alexander, R. C. & Brown, J. C. 2002, Sol. Phys., 210, 407 
Alexander, R. C. & Brown, J. C. 2006, in preparation 
Bai, T. & Ramaty, R. 1978, ApJ, 219, 705 

Bogovalov, S. V, Kotov, Y. D., Zenchenko, V. M., et al. 1985, Soviet Astronomy 

Letters, 11, 322 
Brown, J. C. 1971, Sol. Phys., 18, 489 
Brown, J. C. 1972, Sol. Phys., 26, 441 

Brown, J. C, Emslie, A. G, & Kontar, E. P. 2003, ApJL, 595, LI 15 
Catalano, C. P. & van Allen, J. A. 1973, ApJ, 185, 335 

Conway, A. J., Brown, J. C, Eves, B. A. C, & Kontar, E. 2003, A&A, 407, 725 
Datlowe, D. W., Elcan, M. J., & Hudson, H. S. 1974, Sol. Phys., 39, 155 
Datlowe, D. W., O'Dell, S. L., Peterson, L. E., & Elcan, M. J. 1977, ApJ, 212, 
561 

Drake, J. F. 1971, Sol. Phys., 16, 152 

Elwert, G. & Haug, E. 1971, Sol. Phys., 20, 413 

Farnlk, E, Hudson, H., & Watanabe, T. 1997, A&A, 320, 620 

Johns, C. M. & Lin, R. P. 1992, Sol. Phys., 137, 121 

Kane, S. R. 1974, in IAU Symp. 57: Coronal Disturbances, ed. G. A. Newkirk, 
105 

Kasparova, J., Karlicky, M., Kontar, E. P., Schwartz, R. A., & Dennis, B. R. 

2005, Sol. Phys., 232, 63 
Kontar, E. P. & Brown, J. C. 2006a, Advances in Space Research, 38, 945, 

http://arxiv.org/astro-ph/05084 1 8 
Kontar, E. P. & Brown, J. C. 2006b, ApJL, 653, L149, http://arxiv.org/astro- 

ph/0611170 



Kontar, E. P., MacKinnon, A. L., Schwartz, R. A., & Brown, J. C. 2006, A&A, 
446, 1157 

Kontar, E. P., Piana, M., Massone, A. M., Emslie, A. G., & Brown, J. C. 2004, 

Sol. Phys., 225, 293 
Langer, S. H. & Petrosian, V. 1977, ApJ, 215, 666 
Leach, J. & Petrosian, V. 1983, ApJ, 269, 715 

Lin, R. P., Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Phys., 210, 3 

Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 

McTieman, J. M. & Petrosian, V. 1991, ApJ, 379, 381 

Nitta, N., Dennis, B. R., & Kiplinger, A. L. 1990, ApJ, 353, 313 

Ohki, K.-I. 1969, Sol. Phys., 7, 260 

Petrosian, V. 1973, ApJ, 186, 291 

Phillips, K. J. H. 1973, The Observatory, 93, 17 

Pinter, S. 1969, Sol. Phys., 8, 142 

Press, W., Teukolsky, S., Vetterling, W., & Flannery, B. 1992, Numerical recepies 
in C: the art of scientific computing, 2nd edn. (Cambridge University Press) 
Roy, J.-R. & Datlowe, D. W. 1975, Sol. Phys., 40, 165 

Schwartz, R. A., Csillaghy, A., Tolbert, A. K, et al. 2002, Sol. Phys., 210, 165 
Scullion, E., Kontar, E. P., & Kasparova, J. in preparation 
Smith, D. M., Lin, R. P., Turin, P., et al. 2002, Sol. Phys., 210, 33 
Tomblin, F. F. 1972, ApJ, 171, 377 

Vestrand, W. T., Forrest, D. J., Chupp, E. L., Rieger, E., & Share, G. H. 1987, 

ApJ, 322, 1010 
Zhang, J. & Huang, G. L. 2003, ApJ, 592, L49 
Zhang, J. & Huang, G. L. 2004, Sol. Phys., 219, 135 



